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Abstract 

In this work, we calculate the branching ratios (BRs) and the polarization fractions of sixty two charm- 
less two-body Be meson decays into final states involving one vector and one axial-vector meson (VA) 
or two axial-vector mesons(AA) within the framework of perturbative QCD approach systematically, 
where A is either a ■^Pi or ^Pi axial-vector meson. All considered decay channels can only occur through 
the annihilation topologies in the standard model. Based on the perturbative calculations and phe- 
nomenological analysis, we find the following results: (i) the CP-averaged BRs of the considered sixty 
two Be decays are in the range of 10~^ to 10~®; (ii) since the behavior for ^Pi meson is much different 
from that of "^Pi meson, the BRs of Be — t- A{^ Pi){V, A{^ Pi)) decays are generally larger than that of 
Be — > A{'^ Pi){V, A{'^ Pi)) decays in the perturbative QCD approach; (iii) many considered decays modes, 
such as Be — ^ ai(1260)+a;, 6i(1235)p, etc, have sizable BRs within the reach of the LHCb experiments; 

(iv) the longitudinal polarization fractions of most considered decays are large and play the dominant role; 

(v) the perturbative QCD predictions for several decays involving mixtures of '^Pi and/or ^Pi mesons 
are highly sensitive to the values of the mixing angles, which will be tested by the ongoing LHC and 
forthcoming Super-B experiments; (vi) the CP-violating asymmetries of these considered Be decays are 
absent in the standard model because only one type tree operator is involved. 
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I. INTRODUCTION 



Unlike the Bq meson with q = {u,d,s), the Be meson is the only heavy meson embracing 
two heavy quarks b and c simultaneously. Researchers believe that the Be physics must be 
very rich if the statistics reaches high level. With the running of Large Hadron Collider(LHC) 
experiments, a great number of Be meson events, about 10% of the total B meson data, will 
be collected and this will provide a new platform for both theorists and experimentalists to 
study the perturbative and nonperturbative QCD dynamics, final state interactions, even the 
new physics scenarios beyond the standard model(SM) 

Very recently, we studied the two-body charmless hadronic decays Be — t- PP, PV/VP, VV 
and Be — )• ^Pfhere P, V and A stand for the light pseudo-scalar, vector and axial- vector mesons 
respectively) All these decays can only occur via the annihilation type diagrams in the SM. 

Although the contributions induced by annihilation diagrams are suppressed in the decays of 
ordinary light Bg [q = u, d, s) mesons, they could be large and detected at LHC experiments 0] 
in Be meson decays. According to the discussions as given in Ref. the charmless hadronic Be 
decays with decay rates at the level of 10~^ could be measured at the LHC experiments with the 
accuracy required for the phenomenological analysis, it is therefore believed that they can help 
the people to understand the annihilation decay mechanism in B physics well. In this paper, 
we will extend our previous studies of two-body Be decays to Be — )■ VA and AA modes, which 
are also pure annihilation type decays, and expected to have rich physics since there are three 
polarization states involved in these decays. 

In this paper, we will calculate the CP-averaged branching ratios (BRs) and polarization 
fractions of the sixty two charmless hadronic Be — )• VA, A A decays by employing the low energy 
effective Hamiltonian 0] and the perturbative QCD (pQCD) factorization approach In 
the pQCD approach, the annihilation type diagrams can be calculated analytically, as have 
been done for example in Refs. @, 0,0-14 1. First of all, the size of annihilation contributions 
is an important issue in the B meson physics, and has been studied extensively, for example, 
in Refs. [3, 0, 0, 15,I3]- Second, the internal structure of the axial- vector mesons has been 



one of the hot topics in recent years 174l9l|. Although many efforts on both theoretical and 



experimental sides have been made [2CH.2.6I] to explore it through the studies for the relevant decay 
rates, CP-violating asymmetries, polarization fractions and form factors, etc., we currently still 
know little about the nature of the axial-vector mesons. Furthermore, through the polarization 
studies in the considered Be — >■ VA, A A decays, these channels can shed light on the underlying 
helicity structure of the decay mechanism. 

The paper is organized as follows. In Sec. HH we present the formalism and give the essential 
input quantities, including the operator basis and the mixing angles between ^Pi and/or ^Pi 
mesons. The wave functions and distribution amplitudes for Be and light vector and axial- vector 
mesons are also given here. Then we perform the perturbative calculations for considered decay 
channels in Sec. IIIII The analytic expressions of the decay amplitudes for all considered sixty 
two Be — )• VA, AA decay modes are also collected in this section. The numerical results and 
phenomenological analysis are given in Sec. llVi The main conclusions and a short summary are 
presented in the last section. 



II. INPUT QUANTITIES AND FORMALISM 



In the following we shall briefly discuss the mixing of axial-vector mesons and summarize 
all the input quantities relevant to the present work, such as the operator basis, mixing angles. 
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wave functions and light-cone distribution amplitudes for light vector and axial- vector mesons. 
Finally, the formalism of pQCD approach will also be presented briefly. 



A. Effective Hamiltonian 

For those considered charmless hadronic Be decays, the related weak effective Hamiltonian 
Hcs is given by [H] 

i/eff = ^ [V;,VuD (Cl(/i)Ol(/i) + C2(m)02(m))] , (1) 

with the local four-quark tree operators Oi^2 

O2 = Upl^il - 75)^/3C«7''(l - I5)ba , (2) 

where Vcb,VuD are the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements, "Z?" denotes the 
light down quark d or s and Ci{fi){i = 1, 2) are Wilson coefficients at the renormalization scale 
fi. For the Wilson coefficients Ci{fi), we will also use the leading order expressions, although 
the next-to-leading order calculations already exist in the literature jB]. This is the consistent 
way to cancel the explicit fj, dependence in the theoretical formulas. For the renormalization 
group evolution of the Wilson coefficients from higher scale to lower scale, we use the formulas 
as given in Ref. 0] directly. 



B. Mixtures and Mixing Angles 

In the quark model, there exist two distinct types of light parity-even axial-vector mesons, 
namely, ^Pi ( J^^ = 1++) and ^Pi{J^^ = The ^Pi nonet consists of ai(1260), /i(1285), 

/i(1420) and Kia states, while the ^Pi nonet has 6i(1235), /ii(1170), /ii(1380) and Kib states. 
In the SU(3) flavor limit, these mesons can not mix with each other. Because the s quark is 
heavier than u,d quarks, the physical mass eigenstates i<'i(1270) and iCi(1400) are not purely 
l^Pi or l^Pi states, but believed to be mixtures of Kia and Kib^- Analogous to rj and 77' system, 
the flavor-singlet and flavor-octet axial- vector meson can also mix with each other. It is worth 
mentioning that the mixing angles can be determined by the relevant data, but unfortunately, 
there is no enough data now for these mesons which leaves the mixing angles basically free 
parameters. 

The physical states iri(1270) and 7^1(1400) can be written as the mixtures of the Kia and 
KiB states: 

/ Ki(1270) \^( sinOK cos Ok \ ( Kia\ 

V Ki(1400) J V cos^i^ -sinOK J \ Kib J ^ ' 

The mixing angle 9k still not be well determined because of the poor experimental data. In this 
paper, for simplicity, we will adopt two reference values as those used in Ref. [19I]: 6k = ±45°. 



^ For the sake of simplicity, we will adopt the forms ai, bi, K' , K , f , f , h' and h to denote the axial- vector 
mesons ai(1260), 6i(1235), ifi(1270), i«:i(1400), /i(1285), /i(1420), /ii(1170) and /ii(1380) correspondingly 
in the following sections, unless otherwise stated. We will also use Ki, /i and hi to denote _K'i(1270) and 
-/^i(1400), /i(1285) and /i(1420), and /ii(1170) and /ii(1380) for convenience unless explicitly otherwise stated. 



3 



Analogous to the rj-rj' mixing in the pseudoscalar sector, the /ii(1170) and /ii(1380) (l^Pi 
states) system can be mixed in terms of the pure singlet \hi) and octet \hs), 



/ii(1170) \ _ / cos6li sin6li \ ( hi 
/ii(1380) J ~ V-sin6li cos6li j V /is 

Likewise, /i(1285) and /i(1420) (the l^Pi states) will mix in the form of 

/i(1285) \ / cos 03 sin 03 A ( fi 
/i(1420) J {-sines cos 03 J \ fs 

where the component of hi, fi and /ig, /§ can be written as 



(4) 



(5) 



\hi) = \fi) = -^m + lss)), 

\hs) = \f8) = -^m-2\ss)), (6) 



where q = {u,d). The values of the mixing angles for l^Pi and l^Pi states are chosen as [19l |: 

01 = 10° or 45°; 03 = 38° or 50°. (7) 

C. Wave Functions and Distribution Amplitudes 

In order to calculate the decay amplitude, we should choose the proper wave functions for 
the heavy Be, and light vector and axial-vector mesons. For the wave function of Be meson, we 
adopt the form (see Ref. 0], and references therein) as, 

where the distribution amplitude (pBc would be close to 5{x — mc/mB^) in the non-relativistic 
limit because of the fact that Be meson embraces two heavy quarks. We therefore adopt the 



non-relativistic approximation form for (pBc ^ [27l. l28l| 



(I^bSx) = ^-^=(^(x - mc/msj , (9) 

where Jb^ and Nc are the decay constant of Be meson and the color number, respectively. 

For the wave functions of vector and axial-vector mesons, one longitudinal(L) and two 
transverse(r) polarizations are involved, and can be written as, 

^.(^(x) = -l={mvrvMx)+rvrMx)+rnvct>U^)]e,p , (10) 

$^(x) = -^{mvifct>l{x)+^*^rcf^{x) + mvie,upal^^^e*T^nP^^ , (11) 

c^i(x) = -^j^{mArA<^A{x)+rAr4>A{x) + ^A(t>A{^)],^p , (12) 

$5(x) = -^j,{mArA<Plix)+rAr'pAix) + mAie,.p.75l^^^^^^ , (13) 
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where ^y^^-^ denotes the longitudinal and transverse polarization vectors of vector (axial- vector) 
meson, satisfying P • e = in each polarization, x denotes the momentum fraction carried by 
quark in the meson, and n = (1, 0, Ot) and v = (0, 1, Ot) are dimensionless light-like unit vectors. 
We here adopt the convention e^izs _ -j^ £qj, ^j^^ Levi-Civita tensor e^^'^P , 

The twist-2 distribution amplitudes for the longitudinally and tranversely polarized vector 
meson can be parameterized as: 



(t)v{x) 



3/^ 



:X(1 — X) 



x(l — x) 



1 + Sa'V (2x - 1) + a^^y -(5(2x - 1)^ - 1) 



1 + 3aiy {2x - 1) + aiy -(5(2x - 1)^ - 1) 



(14) 
(15) 



Here fy and fy are the decay constants of the vector meson with longitudinal and tranverse 
polarization, respectively. 



The Gegenbauer moments have been studied extensively in the literatures 29|, |30|] , here we 
adopt the following values from the recent updates 31-13]: 



all^, = 0.03 ± 0.02, a^^, =0.11 ±0.09, a, 
a^K, = 0.04 ± 0.03, a^^. = 0.10 ±0.08, a; 



2p ~ '^2uj 



0.15 ± 0.07, a^^ = 0.18 ±0.08 ; (16) 
0.14 ± 0.06, a^^ = 0.14 ±0.07 . (17) 

The asymptotic forms of the twist-3 distribution amplitudes (j)y and cpy^ are 10|: 



2p — '^2uj 



cP'yix) 



2x-lY, (l)\,{x) 
(1 ± {2x - 1)2), 0» (x) 



(2x-l), 



2V2iV: 
A^/2N, 



{2x-l). 



(18) 
(19) 



The twist-2 distribution amplitudes for the longitudinally and trasversely polarized axial- 
vector and ^Pi mesons can be parameterized as 1^, 24]: 



(t)A{x) 
<fA{x) 



3/ 

^/2iv; 

3/ 



V2iv; 



:X(1 — X) 

■.x{l — x) 



4^ + 3a\^ {2x - 1) + al^ -(5(2x - 1)^ - 1) 
a^j, + 3ai^ {2x - 1) + a^^ ^(5(2x - 1)^ - 1) 



(20) 
(21) 



Here, the definition of these distribution amplitudes (t)A{x) and (l>\{x) satisfy the following 
relations: 



Ct)ip^ {x) 



Pi 



2^/2N\, 



ip^{x} = a, 



Pi 



2y/2Nl 



(22) 



where Opsp^ = 1 and a^ip^ = 1 have been used. 



As for twist-3 distribution amplitudes for axial- vector meson, we use the following form |24i ]: 



<t^\{x) 
cP%ix) 



3/ 



2./2N, 

3/ d 
2y/2Nldx 



aW^x - 1)' ± 2 «i4 (2x - l)(3(2x -lf-l)\, 



[x{l-x){4A + aiA{2x-l))] 



(23) 
(24) 
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3/ 



1 



3/ d 



-x)(a5^ + a)^(2x- 1))} . 



(25) 
(26) 



where / is the decay constant. It should be noted that in the above distribution amphtudes of 
strange axial- vector mesons Kia and Kib, x stands for the momentum fraction carrying by the 
s quark. 

The Gegenbauer moments have been studied extensively in the literatures (see Ref. [l^ and 
references therein), here we adopt the following values: 



''2K, 



-0.02 ±0.02 
-0.04 ±0.03 
-0.07 ±0.04 
0.00 ± 0.26; 
-1.08 ±0.48; 
0.02 ±0.10; 



a 



2K, 



a, 



-1.04 ±0.34 
-1.06 ±0.36 
-1.11 ±0.31 
-0.05 ±0.03 
0.14 ±0.15; 
0.17 ±0.22. 



1.95 ±0.35 
2.00 ± 0.35 
1.95 ±0.35 



4k,a = 0-08 ± 0.09; 



-1.95 ±0.45; 



(27) 



D. Formalism of pQCD approach 

Since the h quark is rather heavy, we work in the frame with the Be meson at rest, i.e., 
with the Be meson momentum Pi = (mBc/\/2)(l, 1,0t) in the light-cone coordinates. For the 
non-leptonic charmless B^ — )• M2M3^ decays, assuming that the M2 (M3) meson moves in the 
plus (minus) z direction carrying the momentum P2 (P3) and the polarization vector €2 (es). 
Then the two final state meson momenta can be written as 



P2 = ^(1 - rl rl Or), Ps = ^(4, 1 - r|, 0^), 



(28) 



respectively, where r2 = ra2/mB^, = m-^/mB^ with 7712 = raM2 and ms = ruM^- The 
longitudinal polarization vectors, €2 and £3 , can be given by 



e'2 



\/2m2 



{l-rl-rlQT), 



mB, 
\/2m3 



{-rll-rlOr). 



(29) 



And the transverse ones are parameterized as e^" = (0,0, It), and = (0,0, It). Putting the 
(light) quark momenta in B^ M2 and M3 mesons as ki, k2, and k^, respectively, we can choose 



A:i = (xiP+,0,kiT), A:2 = (x2P2^,0,k2T), A:3 = (0, X3P3-, k3T). 



(30) 



Then, for Be — >• M2M3 decays, the integration over /c^ , ^2 , and k^ will conceptually lead to 
the decay amplitudes in the pQCD approach. 



^(Pc M2M3) ~ jdxidx2dx3bidbib2db2b3db3 

C{t)^BAxiM)^M2{x2M)^MAxzM)H{xi,bi^t)St{xi)e-^'^^^ 



•Tr 



.(31) 



^ For the sake of simplicity, in the following, we will use M2 and M-j, to denote the final state mesons respectively, 
unless otherwise stated. 
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FIG. 1: Typical Feynman diagrams for charmless hadronic Be — > VA, AA decays. 



where bi is the conjugate space coordinate of kix, and t is the largest energy scale in function 
H{xi,bi,t). The large logarithms ln{mw/t) are included in the Wilson coefficients C{t). The 
large double logarithms (In^Xj) are summed by the threshold resummation 34 1, and they lead 
to St{xi) which smears the endpoint singularities on Xj. The last term, e~^^^\ is the Sudakov 
form factor which suppresses the soft dynamics effectively [35]. Thus it makes the perturbative 
calculation of the hard part H applicable at intermediate scale, i.e., ruB^ scale. We will calculate 
analytically the function H{xi,bi,t) for the considered decays at leading order in expansion 
and give the convoluted amplitudes in next section. 



III. PERTURBATIVE CALCULATIONS IN PQCD APPROACH 

There are three kinds of polarizations of a vector or axial-vector meson, namely, longitudinal 
(L), normal {N), and transverse (T). Similar to the pure annihilation type Be — )■ VV decays 0], 
the amplitudes for a Be meson decaying into one vector and one axial- vector meson or two axial- 
vector mesons are also characterized by the polarization states of these vector and axial-vector 
mesons. In terms of helicities, the decay amplitudes M^'^^ for Be — )■ Af2(-P2; £2)^-^3 (-^3) ^3) decays 
can be generally described by 



a 9^" + 



17121713 



ml Ml + mlMNelia = T) ■ e^(a = T) 
+iMTe''^^PeU^)elp{^)P2^P3p , 



(32) 



where the superscript a denotes the helicity states of one vector and one axial-vector meson or 
two axial-vector mesons with L(T) standing for the longitudinal (transverse) component. And 
the definitions of the amplitudes Mi{i = L,N,T) in terms of the Lorentz- invariant amplitudes 
a, b and c are 



m 



Be 



Ml = a el{L) ■ 4{L) + 



m2m3 



el{L) . P^ el{L) . P2 



171%^ M 



N 

Mt 



(33) 



r2 rs 



We therefore will evaluate the helicity amplitudes Ml,Mn,Mt based on the pQCD factoriza- 
tion approach, respectively. 

In the following we will present analytically the factorization formulas for sixty two charmless 
hadronic Be — )• AV/VA, AA decays. From the effective Hamiltonian ([T]), there are four types of 
diagrams contributing to these considered decays as illustrated in Fig. [1] with single {V — A){V — 
A) currents. From the first two diagrams (a) and (b) in Fig.[Tl by perturbative QCD calculations. 
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we can obtain the Feynman decay amplitudes for factorizable annihilation contributions for 
Be AV, VA, AA as the following sequence, 

rl roo 

Fk{AV) = -STrCpml^ / dx2dx3 / 62^6263^63 

Jo Jo 

X { [x2Mx2)(pvix3) - 2r^r^ ((^2 + 1)0^(2^2) + {x2 - l)^Aix2)) 

XCplr^Xs)] Efa(ta)hfa{l - ^3, X2, 63, 62) + Efa{tb)hfaix2,l - Xs,b2,bs) 

X [(X3 - 1)4>a{x2)4>v{x3) - 24r^<t>'Ax2) ((X3 - 2)0t^(x3) - X3</)V(X3))] } , (34) 

Ff^{AV) = -STTCpml^ / dx2dx3 / 62^6263^63 r^r^ 

Jo Jo 

X {hfa{l - X3, X2, 63, b2)Efa{ta) [{X2 + I)((^i(x2)(^v(x3) + (3^2 ) (a^S ) ) 
+ {X2 - l){rA{x2)rv{x3) + rA{x2Wv{X3))] 

+ [{X3 - 2){4>\{X2)MXZ) + (t>A{x2)Mx^)) - X3 (<^i(x2)</'^(x3) + ^\{x2Wy{x3y)\ 
xEfa{tb)hfaix2, 1 - 0:3, 62, 63)} (35) 

pi POO 

fJ^{AV) = -leTrCpml^ / dx2dx3 / 62^6263^63 r^r^ 

Jo Jo 

X {hfail - X3,X2, 63, b2)Efa{ta) [{x2 + 1) (<^^ (X2)</'^ (0:3) + (X2)</)^ (^3)) 
+ {X2 - l){rAix2mix3) + 4>Aix2Wv{x3))] 

+ [{X3 - 2){rAx2Wv{x3) + rA{X2)rv{X3)) ' X3 (</)i(x2)0?.(x3) + <^^(X2) 0^(0:3))] 
xEfaih) hfa{x2, 1 - X3, 62, 63)} (36) 



/*! poo 

Fk{VA) = -STrCpml^ / dx2dx3 / 62^6263^63 

X { [.X2(/)y(x2)<^A(a:3) + 2r^r^ ((^2 + l)0v'(x2) + (x2 - l)<b\,{x2)) 

^'pA.iXd)] Efa(ta)hfail - .X3 , .X2 , 63 , 62) + E fa{tb)h fa{x2, 1 " ^3,62,63) 

X [{X3 - l)4>v{x2)ct>A{xz) + 2rXri4>Ux2) ((^3 - 2)4>\{x3) - X3ct>\{x3))] } , (37) 



/■l poo 

Ff^{VA) = -SnCFml^ / ^^2^x3 / 62^6263^63 r^r^ 

Jo Jo 

X {/l/a(l -X3,X2,63,62)^/a(ia) [(^2 + l)((/)y (x2)(/)i(x3) + #(x2)</'A(a;3)) 
+ (X2 - l)(<^^(X2)<^i(x3) + rv{X2)4>liX3))] 

+ [(X3 - 2)(^^(X2).^^(X3) + ct>\rix2WA{x3)) - X3 (<^?.(x2).^:i(x3) + <^^(X2)<^^(X3))] 

X-E/„(ib)/i/„(x2,l-X3,62,63)} , (38) 

rl poo 

fUVA) = -IQirCpml^ / dx2dx3 / 62^6263^63 r 2^ 

JO Jo 

X {/i/a(l -X3,X2,63,62)^/a(ta) [(X2 + 1) ((/>?. (x2)(/'a(x3) + (X2)(^A (^^s)) 
+ (X2 - 1)(.^?.(X2).^^(X3) + ct>"y{x2WA{x3))] 

+ [(X3 - 2)(0^(x2)</>:i(x3) + <t>v{x2)<t>\{x3)) - X3 (<^^(X2)<^^(X3) + <^^(x2)0!i(x3))] 
XS/„(ib) /1/„(X2,1 -X3,62,63)} , (39) 
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dx2dx3 / b2db2h^db2 
" Jo Jo 

X { [X2MX2)MX3) + ((^2 + i)<l)2i^2) + (^2 " 1)4(2^2)) 

Xlpl^Xa)] Efa(ta)hfa{l - X3, X2, 63, ^2) + Efa{tb)hfaix2, I - X3,b2,b3) 

X [{X3 - l)Mx2)Mx3) + 24rict)l{x2) ((X3 - 2)(/)^(x3) - xa^^xa))] } , (40) 

f'\ /"OO 

F/^(AA) = SvrCpml^ / ^^2^X3 / b2db2b3dbs 

Jo Jo 

X {/i/a(l - X3,X2,63,&2)^/a(ia) [(x2 + 1) (0^ (X2)</'^ (X3) + 0^ (xa)./.^ (X3)) 
+ (X2 - mU^2W3{x3) + <^^(X2)0I(X3))] 

+ [(X3 - 2)(,/)^(x2)0^(x3) + cPl{x2mXs)) - X3 (0^ (xs)^?^ (X3) + 0^ (X2)</.^ (X3))] 
xEfa{tb)hfaix2,l - X3,b2,b3)} , (41) 

/"l /"OO 

fJ„(AA) = 16^C7^m|^ / dx2dx3 / 62^6263(^63 r-M 

JO JO 

X {/l/a(l -X3,X2,63,62)£^/a(ia) [(x2 + 1) (</'2 (^2)</'3 (^s) + 02 (2;2)</'3 (^^s)) 
+ (X2 - 1)(</)^(X2)</.^(X3) +</.^(x2)0lj(x3))] 

+ [(X3 - 2){r2{x2mX3) + ^I{x2)r3{x3)) - X3 (</'i(x2)</)^(x3) + ./>^(X2)0I(X3))] 
x£'/a(tb) /l/a(x2, 1 - X3,62,63)} , (42) 

where the superscripts V and A in the formulas express the types of mesons involved in the 
considered decays, the subscripts fa and na (to be shown below) are the abbreviations of 
factorizable annihilation and nonfactorizable annihilation respectively, and Cp = 4/3 is a color 
factor. Moreover, the terms proportional to '''2{3) ^'-'^ change the results significantly and 
they have been neglected safely because the values of r"^^^-^ are numerically small: rg^g^ < 5%. 
For the function hfa, the scales U, and Efa{t), we use the expressions as given in Appendix B 
of Ref. 0. 

For the nonfactorizable diagrams (c) and (d) in Fig. [H all three meson wave functions are 
involved. The integration of 63 can be performed using 5 function 5(63 — 62), leaving only 
integration of 61 and 62. The corresponding decay amplitudes are 

1 /n /"OO 
M^„(Ay) = —TrCFml^ / dx2dx3 / 6i(i6i 626^62 

-J JO JO 

X {Ena{tc) [{rc - X3 + 1)0A (x2)</'\/ (X3) - r^r^ ((/.^(X2)((3rc + X2 - X3 + 1) 
X(pvix3) - {rc -X2-X3 + l)</>*^/(x3)) + 0^(x2)((rc - X2 - X3 + l)0y(x3) 
+ (rc - X2 +X3 - l)(py{x3)))] h':^^{x2,X3,bi,b2) - h'^a{x2,X3,biM)Ena{td) 

X [{ri, + Tc + X2 - l)(t)A{x2)4>v{x3) - r2r3 ((/)5i(x2)((4rb + + X2 - X3 - 1) 
X(t>v{x3) - {rc + X2 + X3- l)</.*v(x3)) + (t>A{x2){{rc + X2 + X3- l)(t>v{x3) 
-(r, + X2-X3-l)</.t/(x3)))]} , (43) 

mZ{AV) = -^^TrCpml^ [ dx2dx3 H bidb^b2db2 4r^ 
-J JO JO 

X {rc (X2 ) 0^ (X3 ) + (/)^ (X2 ) (/> v (X3 ) ] £;„a (t c) (2^2 , X3 , 61 , 62 ) 

-r,[(/<^(x2)</'^(x3) + .^:4(x2)0^(x3)]i?na(td)/ii(x2,X3,6i,62)} , (44) 
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-irCpTn^B^ I dx2dx3 / bidbib2db2 r2rY 
" Jo Jo 



3 

X {rc [(Aa (2^2)0^(2:3) + (l)Aix2)(t>vix3)] Ena{tc)Ka{x2,X3, 62) 

-r,,[</.^(x2)</)^(x3) + <^:i(x2)<^^(a;3)]£;„a(td)/ii(x2,X3, 61,62)} . (45) 



16v^ 

M^^{VA) = ^TrCpml^ / dx2 ^^3 / 61^6162(^62 

" Jo Jo 



3 



X {Enaitc) [{rc -X3 + l)0y (x2)</)A(a^3) + (<^t.(x2)((3rc + X2 - X3 + 1) 

X(t^A{x3) - {rc -X2-X3 + 1)(/)*^(X3)) + (t)\^{x2){{rc -X2-X3 + 1)0^(3^3) 

+ (r-c - ^2 + X3 - 1)0^ (X3)))] /1^„(X2, X3, 61, 62) - h^^{x2, X3, bi,b2)Ena{td) 

X [{rb + rc + X2- l)(t)v{x2)(t)A{x3) + r^r^ {^y{x2){{^rb + rc + X2 - ^3 - 1) 
X(t^A{x3) - {rc + X2 + X3- 1)(/>5(X3)) + (t)\,{x2){{rc + a;2 + ^3 - l)(t)A{x3) 

-{rc + X2-X3-l)(t>'A{x3)))]} , (46) 

mZ{VA) = -^^TrCpml^ f\x2dx3 /"°° 61^6162^62 

•J JO JO 

X {rc [(/)^(x2)(/>a(2:3) + 0y(2:2)</>!4(2:3)] Ena{tc)Ka{x2,X3, 61, 62) 

-r6[0?^(x2)0i(x3) + </'^(x2)0:i(x3)]i?na(id)/ina(x2,X3,6i,62)} , (47) 

M^^{VA) = TrCpml^ / dx2dx3 / 61^6162^62 r^r^ 

•J Jo Jo 

X {rc [(/)^(X2)(/)^(X3) + #(x2)</)a(2;3)] -Ena(tc)/ina(a;2,a;3,6l,62) 

-rb[<l>v{x2)<j>A{x3) + (Pv{x2WA{x3)]Ena{td)hL{x2,X3,bi,b2)] • (48) 

16\/6 /"^ /"°° 

M^„(^^) = — — 7rCFm|^ / dx2dx3 / 61^6162^62 
^ " Jo Jo 

X {^na(tc) [(rc - X3 + l)Mx2)Mx3) + r^^^ ((/>i(a:2)((3rc + a;2 - X3 + 1) 

X(t)l{x3) - {rc -X2-X3 + l)(j)i{x3)) + (i)i{x2){irc -X2-X3 + l)</'3(-T3) 
+ (rc- X2 + X3 - l)(/)|(a;3)))] h^„^{x2,X3,bi,b2) - h!^aa{x2,X3,biM)Ena{td) 

X [(rfe + rc + X2 - l)(/'2(a;2)</'3(a;3) + r^r^ ((/>2(x2)((4rf, + rc + X2 - X3 - 1) 

X03(2;3) - {rc + X2+ X3 - l)(t>\{x3)) + 4)\{x2){{rc + X2 + X3- l)(t>i{x3) 

-{rc + X2-X3-l)(ti{x3)))]] , (49) 

MZ{AA) = ^^nCFml^ [\x2dx3 H bidbib2db2riri 
J Jo Jo 

X {rc [(l>^{x2)(l>t{X3) + <^^(X2)01(X3)] ^na(tc)^na(a^2,X3,6i,62) 

-rb[r2{x2)4>U-^3) + ^l{x2W3{x3)]Ena{td)ht{x2,X3M,b2)} , (50) 

M^^{AA) = ——TTCFml^ / dx2dx3 / 61^6162^62 r^r^ 
J Jo Jo 

X {rc [05(X2)<^^(X3) + 4>l{x2)4>t{x3)] Ena{tc)Ka{X2, X3MM) 

-r6[</.5(x2)</.l(x3)+</.^(x2)0g(ar3)]-E„„(td)/i^„(x2,X3,6i,62)} . (51) 

where r^, = mh/niB^, rc = rnc/mB^, and r^ + rc ~ 1 for Be meson. 
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There are three kinds of polarizations in these Be — )• VA, AA decays, namely, longitudinal 
(L), normal {N) and transverse (T). The decay amplitudes are classified accordingly, with 
H = L,N,T. Prom the effective Hamiltonian ([T]), based on Eqs. (j34ti5ip . we can combine all 
contributions to these considered decays and obtain the total decay amplitude generally as. 



(52) 



where ai = Ci/3 + C2. Then we can write down the total decay amplitudes for sixty two 
charmless two-body nonleptonic Be meson decays into final states involving one vector and one 
axial-vector meson (VA) or two axial-vector mesons{AA) one by one. 



I. Be —i- VA/AV decay modes 
(i) For A5 = processes, 

V2M^{Bc^ 



V2M" (Be ^ a+ 



V2M"{Be at 00) 



{ 



cb^ud 



v:,v^d ' 



fB.FZ.1;a, + M2i;c, 



fa;H 

fB.FjQa, + M2^Cr 



(53) 



(54) 



(55) 



V2M^ {Be ^ b+ 



V:,V.d { [fB.F0Ha, + M^%C^ 
-V2M^{Be ^ p+b?) 



v:,v^d 

+ 



fBF)]:H<^i + M:^-C, 



fB.F^^.Jjai + M^^.jjCi 



(56) 



(57) 



(58) 



M^'iBe^p+f) = V:,Vud 



cos ^3 



fB.iFfJn + FfAhl 



smyg 
H — [jBc 



iFfafn + Fji^)ai + (M^ff^ + M£f^)Ci] } , 



(59) 
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sm (^3 



V3 

HM^nJ„ + M%)C, 



cosd^a 

7^ [/Be 



V6 



cos 6*1 



V3 



H 7^ [/Be 



V6 



— sm f^i 



x/6 



+ cos^.[/.,FS--a,+M„^--C,]}/ 



+ cos Ok 

M''{Bc^K^K*+) = KfeKdjcos^^ 

— sin 6 k 



fy-0 K* K* 



fB.F^If*<^l+KTH^ Cl 



For AS = 1 processes, 
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+ cos Ok [fB.F^^fa, + ul^ffC,] } , 



V^V^s {cos Ok [fB.F^fa^ + Mfjf Ci | 



sin 



na;H ^1 



}■ 



iVus {cos 9k [fB.F^Ifa^ + uf^-^C, 



cos ^3 



+ J 

VD 



■ sm t)s 



cos 6/3 
H — [jBc 



H 7^ [/Be 



M"{Bc^ K*+h") = 



V*V J_^^ 



fBAFl.!'H^+FYJ^)a, 



COS^l 

H — [JB^ 



a/6 
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M^{B^ <t)K"' 



5c -> 74^4 decay modes 
(i) For AS = processes, 



\/27W^(5c^a+5?) 



M"{B, a+/') 



V6 



</8 
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V3 

Mfi ^ njftf't.r. 1 , sink's 



[fE 



7J1 "1 



)ai 



M"{Bc a+h') 



/ pdi rag 



M^iB^^bth') = V:,V^,[^[fsAF^i;S+F';l$)a, 



h+h" h<ih+ sin ^1 



sin ^1 h+^i" 



+(KQ+<S)Ci] + 



V6 
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+ cos2e;,(/s.Fj||°^-ai + M5f°^-Ci)| , (91) 



-cos0Ksin0;,(/B.Fjj|"^-ai + M5^°^-Ci)| , (92) 



+ sin^ 6K{fB.Fj$;^--a, + M^l^^-Ci) 

+ cose^sin^x(/B.i^J|l"''^"ai + M5f"^-Ci)| , (93) 

>t^(Be^F^if"+) = y;,K,|cos2^^(/5^Fj!^"^i^ai + M5^"^i^Ci) 

-cos0Ksin^^(/B,Fj|^°^-ai + m5^°^-C7i) 
+ cos0Ksin^^(/B.Fj!|"^-ai + m5^°^-Ci) 
-sin20^(/5^Fjj|°^-ai+Mj^°^-Ci)| . (94) 

(ii) For AS = 1 processes, 

+ cos9K[fB.F^Jf^a^ + M;^J-/c^]^ , (95) 



-sin^;,[/B.F;^|r^ai + Mf,?-»^Ci]} , (96) 
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+ cos eKifs.F^Jfa, + Mf^fc] 



f h 



(97) 



(98) 



+ i^na;i? + ^na;H J<^lJ + UB 

fpKiAfs OV^^^'^^'\n J-(M^'^^^^ 9A/f■^8■'^l'*^ 

■'-'lJ + ^/l^ UBA^fa;H +^fa;H )«1 

, .^i^iB/i" , ^/fi^iBx^ , cos^x sin 6*3 . 

+ l^na;if + ^na;H + L/Bc 

•(^/^^'^"-2i^S")ai + (Mf„-^ , (99) 

[ -sin 03 singj^ , K^Af^ fiK^A. 
Vch^us<y [JBA^fa-H + ^ fa;H )«1 

, /,^^iA/r , A^/f^iAN^ 1 , COS 03 sin 0g 

+ K^na;H + ^na;H )<-lJ + ^ UBc 

COS sin 03 ^p^is/i" 

■<-^lJ UBA^fa;H + ^ fa;H )«1 

, .^(.i^iB/r , 1 , COS0i£COS03 
+ (^na;i/ + ^na;H )<^lJ + [JBc 

iFf:f - ^Fli^^ai + {M;:.f - 2mS-)C,]} ,(100) 
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fa;H fa;H ^"1 ^ l-'^-'na;/f ^''^^na;// 

COS 03 sill (?A- (pl<lBfr , j^flK.B^ 
■<^lJ UBA^fa;H + ^ fa;H )«1 

+ (^na;ff + ^na;H )<^lJ y| i/Sc 

•(^a!^'^" - 2i^S^)ai + {M^:.f - 2m£5-)Ci]} ,(101) 



•<^lJ + UBA^fa;H + ^ fa;n )«1 

, / , .i^iB/r , Mf^^^^\r 1 COS ^3 sin r . 

+ (^na;/f + ^na;H i/B= 

•(^S^'^" - 2i^/S^")'^i + i^aTn' - 2mS-)Ci]} ,(102) 



+(Mf„:^^" + M„^I--)C7,] + !^^^^ 

fa;H fa;H ^"1 V^'^na;// '^^^^na;H > 

+ UBA^fa;H + ^ fa;H ^ 

iFfX' - '^F^li^a, + (M„^„;^« - 2M„^f^^-)Ci]} , (103) 
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■<^lJ iJBA-t'fa-H +^fa;H )«1 

, , , ,^fKisx^ 1 , COS^^COS^ 

- 2i^S^")«i + - 2M^f J-)Ci]} , (104) 



x/3 

+ i^na;if + ^na;i? )<-^lJ + iJBc 

■<^lJ UBA^fa;H + ^ fa;H )«1 

■K-f' - ^F^Hinai + (M„^„;^«" - 2M„^f^^-)Ci]} , (105) 



+ \^na;H + ^na;H )<^li + UBc 

sin 03 sin r . , j^k^bK , p'^f -^is ^„ 

•<^lJ + UBA^fa;H + ^ fa;H ^ 

^f/W^i^'^i + M'^i-^i^^rl COS 03 sin 

- 2i^;i5^")«i + - 2M„^|5-)Ci]} .(106) 



IV. NUMERICAL RESULTS AND DISCUSSIONS 

In this section, we will calculate numerically the BRs and polarization fractions for those 
considered sixty two Be — >■ VA/AV, AA decay modes. First of all, the central values of the input 
parameters to be used are given in the following. 
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Masses (GeV): 



mw 


= 80.41; 




= 6.286; 


mb 


= 4.8; 


nic -- 


= 1.5; 




= 1.02; 




= 0.892; 


mp 


= 0.770; 




= 0.782 




= 1.23; 




= 1.32; 




= 1.28; 




= 1.28; 




= 1.21; 




= 1.34; 




= 1.23; 




= 1.37 



(107) 



Decay constants (GeV): 



fcf) 


= 0.231; 


fj — 


0.200; 


fK' = 


-- 0.217; 


/I. 


= 0.185; 


fp 


= 0.209; 


f = 
Jp 


0.165; 


fu-- 


= 0.195; 




= 0.145; 


fai 


= 0.238; 


fKiA 


= 0.250; 


fh 


= 0.245; 


ffs 


= 0.239; 


fbi 


= 0.180; 




= 0.190; 


fhi 


= 0.180; 


fh 


, = 0.190 




= 0.489; 















QCD scale and Be meson lifetime: 



A. 



MS 



0.250 GeV, 



0.46 ps. 



(108) 



(109) 



For the CKM matrix elements, here we adopt the Wolfenstein parametrization for the CKM 
matrix, and take A = 0.814 and A = 0.2257, p = 0.135 and fj = 0.349 [2^. In numerical 
calculations, central values of input parameters will be used implicitly unless otherwise stated. 

For these considered Be — )• M2M3 decays, the decay rate can be written explicitly as, 



G'p\Pc\ 

Be a=LT 



(110) 



|P2z| = iPszl is the momentum of either of the outgoing vector or axial- vector 



where |Pc| 
mesons. 

Based on the helicity amplitudes (j33|) . we can define the transversity amplitudes, 



Be 



M 



N, 



A 



1)Mt 



(111) 



for the longitudinal, parallel, and perpendicular polarizations, respectively, with the normal- 
ization factor ^ = G|,Pc/ (167rm^^r) and the ratio r = P2 ■ P-i/{m2 w.3). These amplitudes 
satisfy the relation, 

|^l|2 + |^I,P + |^x|' = 1 (112) 
following the summation in Eq. piO|) . 



The polarization fractions and /_l can be defined as. 



\Aj 



\Al? + \A\\ 



\A, 



(113) 



With the analytic formulas for the complete decay amplitudes as given explicitly in Eqs. (|53t - 
I106p . by employing the input parameters and Eq. llllOp . we calculate and then present the pQCD 
predictions for the CP-averaged BRs and longitudinal polarization fractions (LPFs) of the con- 
sidered decays with errors in Tables HEl The dominant errors come from the uncertainties 
of charm quark mass rric = 1.5 it 0.15 GeV and the combined Gegenbauer moments at of the 
relevant meson distribution amplitudes, respectively. 

Based on the numerical results as given in Tables HllXl some remarks are in order: 
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Among the considered sixty two pure annihilation Be — )• AV/VA, AA decays, the pQCD 
predictions for the CP-averaged BRs of those AS* = processes are generally much larger 
than those of AS* = 1 channels (one of the two final state mesons is a strange meson) , the 
main reason is the enhancement of the large CKM factor ~ 19 for those A^* = 

decays as expected in general. Maybe there exists no such large differences for certain 
decays, which is just because the enhancement arising from the CKM factor is partially 
cancelled by the difference between the magnitude of individual decay amplitude. 

There is no CP violation for all these sixty two decays within the SM, since there is only 
one kind of tree operator involved in the decay amplitude of all considered Be decays, 
which can be seen directly from Eq. (j52p . 



For the ten Be — s- {ai,bi)V decays, the pQCD predictions of the BRs and LPFs for 
both AS = and AS = 1 processes are listed in the Tabled! As argued in Ref. [3], 
the LHCb experiments could observe the BRs of annihilation Be meson decays at the 
level of 10~^, the decays Be — )■ aiio,bip will thus be detected at LHC because they are 
just within its reach. As for the polarization, all these ten decays are governed by the 
longitudinal contributions. The LPFs are around 95% within the theoretical errors except 
for Be — )• afuj,aiK* modes (~ 85%), and for Be — )• aip channels, fiiBe aip) ~ 1, 
which will be tested in the LHCb experiments. 

Since the behavior of hi meson is contrary to that of ai meson, one can find that Br {Be — )• 
hip) > Br{Be — )■ dip) as given in Table [I] and the ratio of the corresponding BRs for 
Be — )■ aip and Be — )■ bip is that 

BrjBe ^ b+/) ^ Br{Be ^ b\p+) ^ ^ 

Br{Be atp^) Br{Be alp+) ^ ' ' ^ > 

Similarly, for Be — t- aiK* ,biK* decays, the BRs of the latter modes are larger than that 
of the former ones and 

Br(Be ^ h+K*^) ^ BrjBe ^ 6?^*+) ^ ^ 

Br{Be a^K*^) Br{Be a^K*^) ^ ' ' ^ ' 

The above two ratios exhibit the annihilation decay pattern consistent with those as shown 
in Ref. EF 



TABLE I: The pQCD predictions of BRs and LPFs for Be {ai,bi){p, K*,uj) decays. The source of the 
dominant errors is explained in the text. 



AS^O 
Decay modes 


BRs (10-'^) LPFs (%) 


AS^O 
Decay modes 


BRs (10-S) LPFs (%) 


Be -+ p+al 
Be ^ afpO 


7.5t°-lim,)tliia.) 99.21^-3 
7.5t°-i(me)l^;0(a,) 99.2tU 


Be ^ p+b? 
Be ^ b+p" 


9.9lt^(me)ll:^(a.) 93.8^ 


Be -> a^ui 


20.2tg:°(me)trj(a,) 84.71^;° 


Be — > h^LO 


0.6l°i(me)l°i(a,) 96.3l?J 


AS=1 
Decay modes 


BRs (10-8) LPFs (%) 


AS^^l 
Decay modes 


BRs (10-^) LPFs (%) 


Be ^ a+K*° 
Be -> K*+a? 


3.3;°:J(me)li:^(a.) 83.6t?j 


Be ^ b+K*° 
Be ^ K* + b? 


3.6lJ:i(me)li:«(a.) 91.2±2-5 
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TABLE II: Same as TableUbut for Be (ai, 6i)(ai, 61) decays. 



AS' = 






AS = 






Decay modes 


BRs (10-5) 


LPFs (%) 


Decay modes 


BRs (10-5) 


LPFs (%) 




0.0 




Be ^ b+b? 


0.0 




Be ^ a+b? 




92Atli 


Be ^ b+a? 


2.2l°j(™e)tJi(aO 


91.8t|° 



• Analogous to Be — >• p'^'/f'^ decay [2] , the contributions from uu and components cancel 
each other exactly and result in the zero BRs for Be — >■ afa^ and Be — >■ bft^. Any other 
nonzero data for these two channels may indicate the effects of exotic new physics. While 
for Be — )• afb^ and Be ^ b clS expected from the analytic expressions, Eqs. (j81|82p . 
due to the same component of uu — dd involved in both axial- vector a\ and b\ mesons at 
the quark level, the pQCD predictions for the BRs and LPFs as given in Table |TI] show 
the identical results as they should be, 

Br{Be alb\) = Br{Be b"[a\) ^ 2.2 x 10^^ , 
fUBe ^ a+6?) = hiBe ^ 6+a?) « 92% . (116) 

where the large BRs (~ 10"^) are within the reach of the LHCb experiments y] and could 
be detected at LHC. 

• Since the ^Pi meson behaves like the vector meson and fai ~ fp from Eq. pOSp . the 
pQCD predictions of BRs exhibit the good consistency generally for Be — )■ afu and Be — ?■ 
p+uj, Be a+K*^ and Be p+K*^, Be a^6?(a06+) and Be p+fe?(p°fe^) decays, 
respectively, within the theoretical errors as roughly estimated. As for the polarizations, 
which can well manifest the helicity structure for the corresponding modes, the LPFs 
present the different features from the decay rates except for Be — >• afb^{a^bf) and Be — s- 
p+b\{p%'[) decays. From Table H the LPFs for Be —7- a^oj and Be — s- a^K* can be read 
straight forward as: /l(5c ^ a^w) = (84.71^;^)% and /l(Bc ^ a+K*^) = (83.6;^|)%. 
As given in Ref. 0, /l(S, ^ p+w) = (92.9l2;0)% and /l(B, ^ p+K*"^) = (94.9t2;0)%, 
where the various errors as specified have been added in quadrature. The above results 
and discussions would be tested with high precision by the relevant experiments operated 



TABLE III: Same as TableUbut for (p, is:*)(/i(1285), /i(1420)) decays. 



AS = Q 
Decay modes 


6I3 = 38° 
BRs {IQ-^) LPFs (%) 


03 = 50° 
BRs(10-'^) LPFs (%) 


Be ^ p+fi(1285) 

Be p+fi(1420) X 10" 


2.lt;j:i^(me)tnK) 82.6i5-J 
0.4tn(-c)t^-j(«.) 88.7t^i,« 


l.^tlii^cniia.) 82.3t5-6 
2.2tr,{m.^tll{a,) Sb.ltH 


A5 = 1 
Decay modes 


6I3 = 38° 
BRs (lO-'^) LPFs (%) 


63 = 50° 
BRs (10-7) LPFs (%) 


Be ^ K*+fi(1285) X 10 
Be ^ K*+fi(1420) 


i.6t;j:i^(me)t;,j(a,) ei.oi^^J 

l.ltEJ:i(me)lHK) 85.4t4.J 


0.4l°:J(™e)t^J(aO 33.7l:lli 
1.2tO;2(me)t;^i{«i) 83.9^5-2 



"Here, the factor 10 is specifically used for the BRs. The following ones have the same meaning. 
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at the ongoing LHC and forthcoming Super-B to identify the hehcity structure even decay 
mechanism in these considered channels. 



In Table IIIH from the pQCD predictions of the BRs and LPFs for Be — )• 
/9+ (/i (1285), /i (1420)) (A5 = 0) and B, ^ ir*+(/i(1285), /i (1420))(A5 = 1) decays, 
one can observe that the BRs of Be — >■ /9"*'/i(1420), -ftr*^/i(1285) are more sensitive than 
those of Be p+/i(1285), K*+/i(1420) to the mixing angle 6*3, 



Br{Be^ p+fiiU20))\ 



6*3=50° 



Br{Be^ p+h{U20))\es=38° 
Br{Be^K*+h{1285))\e,=3S^ 



Br{Be^ K*+fi{1285))\ 



6*3=50° 



5.5 ; (117) 
4.0 ; (118) 



These two relations, Eqs. (I117|118p . can be understood as that the interferences between 
Be —!- p+/i(-Bc K*'^ fi) and Be p+fsiBe K*^fs) become highly construc- 
tive(destructive) to Be p+fi{U20){Be K*+/i(1285)) with the mixing angle 6*3 
changing from 38° to 50°. Moreover, 

Br{Be ^ p+ fi{l285)) [52.5 for 03 = 38° 



Br{Be ^ p+ fi{U20)) [8.6 for ^3 = 50° 

■Br(^c^i^*+/i(1420)) _ f 6.9 for 6^ = 38° 
Br{Be ^ K*+ fi{1285)) ^ 1 30.0 for 6*3 = 50° 



(119) 
(120) 



From the decay amplitudes as given in Eqs. (|59l60l73|74p . the above two relations 
can be understood as follows: (a) for Be — )• /o"^ (/i (1285), /i (1420)) decays, the mix- 
ing coefficients for the former decay are cos03 and sin03, while that for the latter one 
are — sin03 and cos03. For the common component qq, it is found that the contri- 
butions from Be —7- p~^ fi and Be — s- p^ fs interfere constructively(destructively) for 
Be p+fi{128b){Be p+/i(1420)); (b) for Be i^*+(/i(1285), /i(1420)) channels, the 
mixing parameters remain unchanged, however, a new part of contribution from ss compo- 
nent involved in both fi and /§ with different signs results in the construction (destruction) 
to Be K*+fi{U20){Be K*+/i(1285)). Additionally, the LPFs for these decays are 
stable to the mixing angle and play the dominant role except for fiiBe — s- i^*'^/i(1285)), 
whose value change from 61.0% at ^3 = 38° to 33.7% at ^3 = 50°, which will be confronted 
with the relevant experiments in the future. 



TABLE IV: Same as Tabic Ubut for {p, K*){hi{1170), hi{1380)) decays. 



A5' = 
Decay modes 


01 ^ 10° 
BRs (10-^) LPFs (%) 


61 = 45° 
BRs (10-^) LPFs (%) 


Be -+p+hi(1170) 
Be ^ p+hi(1380) 


6Atti{m,)tli{a.) 96.411;^ 
2.5tli{m,)ll"M) 96.3t2;i 


0.3t°J(me)t°:?(a.) 98.0t^] 


AS'= 1 
Decay modes 


01 ^ 10° 
BRs (10-^) LPFs (%) 


61 = 45° 
BRs (10-^) LPFs (%) 


Be ^ K*+hi(1170) 
Be ^ K*+hi(1380) 


0.2ig:i(me)inK) 82.3t^-j 
4.8t?-j(me)l^:^o(a.) 93.2^11 


2.lt°:^(m,)tJi(a,) 89.5t^-4 
2.9t\l{m,)t\\{a,) 9^.9t\l 
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• The pQCD predictions for Be — )• {p'^,K*~^)hi decays, as given in Table HV] can be ex- 
plained in a similar way as for Be — )• {p'^ , K*~^) fi. 

• The numerical pQCD results for AS' = 5^ ^ (a]^, 6+)(/i(1285), /i(1420)) decays, 
as given in Table |Vl can be commented in order: (a) the BRs of these modes depend 
weakly on the mixing angle ^3 except for Be — )• a]''/i(1420); (b) these 4 considered de- 
cays are governed by the longitudinal contributions for both 6*3 = 38° and ^3 = 50°; 
(c) as mentioned in the text above, l^Pi meson behaves close to the vector meson, the 
phenomenology of B^ — >• a^(/i(1285), /i(1420)) can therefore be understood as that of 
Be — )• (/i (1285), /i (1420)); (d) the mixing factors cos 9^ and sin 6*3 make the inter- 
ference between Be — )• {af,bf)fi and Be — )■ {a^,b^)fs constructive (destructive) to 
Be ^ (a+,6+)/i(1285)(i?e ^ (a+, 6+)/i(1420)). 

• From the predictions for Be — )• (a^, 6^)(/ii(1170), /ii(1380)) presented in Table IVTl 
some discussions could be addressed as: (a) since the behavior of l^Pi meson is dif- 
ferent even contrary to that of l^Pi meson, a surprisingly large branching ratio for 
Be — ^ 6^/ii(1700)(~ lO"'^) with the constructive effects induced by the interference be- 
tween Be — )• b^hi and Be — )• b^hs is produced, which will be tested stringently by the 
forthcoming relevant LHC experiments; (b) once the large BRs are verified by the measure- 
ments, the mixing angle Oi could be well determined to improve the precision of the pertur- 
bative calculations; (c) except for Be — )• bf hi{1170) decay, the rest three channels are sen- 
sitive significantly to the mixing ang le 01] (d) similar to Be (oj*", 6|)(/i(1285), /i(1420)) 
decays, the longitudinal components play the dominant role for these 4 channels. 

• For the AS = Be^ K* and Be K^K*^ decays, one can see from Table IVTIl that 
the BRs are large in the range of 10~^ ~ 10~^, which can be detected at the ongoing LHC 
and forthcoming Super B experiments. Moreover, the corresponding ratios of the BRs for 
these considered channels 



Br{Be K*^ Ki{l270)+) Br{Be Ki{l270fK*- 



Br{Be ^ ;^Ki(1400)+) Br {Be ^ i^i(1400)0K*+) 
f 01 9k = 45°, while 

Br{Be^K*°Ki{1270)+) _ Br{Be^Kl{l270fK*^) ^ 1 

Br{Be^K^Ki{Um)+) ~ Br{Be ^ :^(1400)0K*+) ^ L7 



1.7, (121) 



(122) 



TABLE V: Same as TableUbut for (a+, 6+)(/i(1285), /i(1420)) decays. 



A5 - 
Decay modes 


6I3 = 38° 
BRs (10-S) LPFs (%) 


03 - 50° 
BRs (10-^) LPFs (%) 


Be ^ ai(1260)+fi(1285) 

Be ^ ai(1260)+fi(1420) x 10 


6.5tJ-^(me)t;-g(aO 83.61^-4 
Q.itrMc)tll{a.) 56.8t^^j 


6.i±;,:^(me)l°:^9(a.) 84.o;^:g 

3.9tS:^(me)li-J(a,) TS.Stlg^ 


A5 
Decay modes 


6I3 = 38° 
BRs (10-^) LPFs (%) 


6k = 50° 
BRs (lO-'^) LPFs (%) 


Be ^ bi(1235)+fi(1285) 
Be ^ bi(1235)+fi(1420) 


2.8t^i(me)tJ:«(a,) 65.2™ 
1.4t°:2(me)t;j:^(a,) 100.0 ±0.0 


3.0t^:^(me)lJ:^(a,) Q8.7tl\l 
1.2;°;l(me)tJi(a,) 100.0t°:° 
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TABLE VI: Same as TableUbut for (a+, 6+)(ft.i(1170), /ii(1380)) decays. 



A5' = 
Decay modes 


9i = 10° 
BRs (10-*^) LPFs (%) 


6*1 = 45° 
BRs (10-^) LPFs (%) 


Be ^ ai(1260)+hi(1170) 

Be ^ ai(1260)+hi(1380) x 10 


1.3±°:^(me)t°i(a,) 86.3i^:° 
l.ltroirnc)tliia.) 68.8^1^ 


0.7lg:l(me)l°:^(a,) 73.1^^3^3 
6.8t?:?(me)l^i(a,) 100.0t?;° 


AS^O 
Decay modes 


01 = 10° 
BRs (10-5) LPFs (%) 


6li = 45° 
BRs (10-5) LPFs (%) 


Be ^ bi(1235)+hi(1170) 
Be ^ bi(1235)+hi(1380) 


8.ltliim,)tli{a,) 96.4t}:° 
2.5t°-?(™e)t}1(a0 lOO.Oin 


10.3tti{m,)tt-lia^) 96.4t?i 
0.3tJ] j(me)lo|(aO 100.0 ±0.0 



for Ok = —45°, which indicate that one could determine the size and sign of the mix- 
ing angle 6k after enough Be events become available at the LHC experiments and then 
improve the precision of the theoretical predictions. In terms of polarization, the lon- 
gitudinal contributions play the dominated role for both 9k = 45° and Ok = —45° in 
Be 'K^K*+ modes. In the B^ ^(i^i(1270)+, i^i(1400)+) decays, the transverse 
components govern the former channel for 6k = —45° while dominate the latter one for 
Ok = 45°. These results will be tested by the relevant measurements in the future. 

• Form the numerical results for Be — s- K^{p, lo), the AS = 1 processes, as displayed in the 



TABLE VII: Same as TableUbut for ^ (i^i (1270), i^i (1400)) (p, i^T*, decays. 



AS^O 
Decay modes 


Ok = 45° 
BRs (10-S) LPFs (%) 


6k - -45° 
BRs (10-6) LPFs (%) 


Be ^K*''Ki(1270)+ 
Be ^ K*°Ki(1400)+ 


3.8l^i(me)l^:^(a.) 96.8lii 
2.2l°i(me)l?:°(a.) 42.71^^:^ 


2.3t°i(me)t?:^(a0 42.3tiS-| 
3.8t°|(me)t^:^(a.) 96.9t^;3 


Be ^ KT(1270)°K* + 
Be ^ KT(1400)°K* + 


9.7tliimc)n-'M) 97.5li-J 
5.8tl3imc)tt'M) 82Al:li, 


5.6t|t(™c)t^j(a0 81.ltl°i 
9.6t|^(me)t^:?(a0 97.6t^:^ 


AS^l 
Decay modes 


6k - 45° 
BRs (10-^) LPFs (%) 


= -45° 

BRs (10-^) LPFs (%) 


Be ^ Ki(1270)V 
Be ^ Ki(1400)V 


3.ll}-t(me)l?-?(a,) 89.515-^ 
4.0i;:^(me)l|^(a.) 99.11^4 


4.0t}:^(me)t|?(a.) 99.lt°:« 
3.0t}:^(me)t?j(a0 89.7t5:6 


Be ^ Ki(1270)+p" 

Be ^ Ki(1400) + /9" 


1.5l°:^(me)lJ:^7K) 89.515-J 
2.Qt',%{m,)t\l{a.) 99.ll°-i 


2.0tli{m,)t\l{a,) 99.lt°J 
1.5t°:I(me)tJ:^(aO 89.7t5;6 


Be ^ Ki(1270)+a; 
Be ^ Ki(1400)+a; 


l-4tl]-i(me)tJ:?(a,) 89.8tl'r 

l-7tSj(me)tg:^(a,) 99.i;°:« 


l-7t^:^(me)t°:^(a.) 99.lt°:« 
1.4i°j(™e)±J:^(aO 89.9t5j 


Be ^ Ki(127O)+0 
Be ^ Ki(1400)+(/) 


1.9tS]-j(me)t}:°(a.) 9^.2tlL 
l-4tS1(me)t^,-j(a.) 30.31^1-0 


l-5t[!i(mc)t?,:I(aO 29.9tf,:« 
1.9t°i(™e)t}:^(a0 95.2t?„% 
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Table IVTH one can straightforwardly observe that 

MB, 

-Br{Bc^ Kiil270fp+); (123) 



Br{Bc^ Ki{1270yuj) ~ Br{Bc ^ Ki{l270)+ p^) 

1 



2 

Br{Bc^ Ki{UOO)+oj) ~ Br{Bc ^ Ki{UOO)+ 

= ^Br{B, ^ Ki(1400)V+) ; (124) 

/L(Se^Ki(1270)+^) ~ /L(5,^Ki(1270)+pO) 

= fL{B,^Ki{1270fp+) ; (125) 

fL{Bc^Kiiim)+uj) - /i(5e^Ki(1400)+p°) 

= fUB,^ Ki{U00fp+) . (126) 

within errors for both 9k = 45° and 6k = —45°, where the longitudinal components 
contribute to these considered decays dominantly. The pattern of these decays shown in 
Eqs. () 123111261) can be understood as follows: only the same component uu in both of p^ 
and LO mesons contributes to these physical observables, where the differences mainly arise 
from the different decay constants. Furthermore, by comparison with Be K* and 
Be — )• KiK*^ decays, one can find that the pQCD predictions of Be — )• K^{p,uj) show 
the weak dependance on the value of the mixing angle 9k, which will also be tested by 
the LHC measurements. 

For Be — >• K^(/) decays, it is interesting to note that the decay rates, as listed in Ta- 
ble VII, are close to each other within the theoretical uncertainties, however, the LPFs 



TABLE VIII: Same as Table Ubut for Be ^ (i^i(1270), iri(1400))(ai, 6i, ifi(1270), ifi(1400)) decays. 



AS = 
Decay niodes 


9k = 45° 
BRs (10-5) LPFs (%) 


9k = -45° 
BRs (10-5) LPFs (%) 


Be ^ KT(1270)OKi(1270)+ 
Be ^ KT(1270)°Ki(1400)+ 


1.2t°;?(m,)ljj(a0 99.7t?:J 
3.7tli{m,)tl-l{a,) 96.2tli 


2.9tl■l{me)tttia^) 71.91^°:^ 
1.9l°:5(me)l?:l(aO 94.8l?o^3 


Be -l>Kr(1400)°Ki(1270)+ 
Be KT(1400)OKi(1400)+ 


2.8tl-',{me)tti{ai} 72.7tlli 


3.7ll:?(me)l^:?(a0 96.11^^ 
l.ltVoime)tl-U-^) 99.71?:° 


AS-^ 1 
Decay modes 


Ok = 45° 
BRs (10-^) LPFs (%) 


9k = -45° 
BRs (10-^) LPFs (%) 


Be ^ Ki(1270)'^ai(1260)+ 
Be ^ Ki(1400)°ai(1260)+ 


4.6t};ij(me)t^;I(aO 79.2t\li 
8.0t};?(TOe)l^;'7(«0 100.01°:° 


S.3t\iim,)tl:Ua^) 99.3l°:« 
4.5l}:?(me)l^:^(a.) 81.311^:5 


Be ^ Ki(1270)+ai(1260)0 
Be ^ Ki(1400)+ai(1260)° 


2.31q 5(mc)l]^ 3(0^) 79.2l-^g 3 
^.Ot'oi{m.)tl:Ua.) 100.01°:° 


4.2l?:^(™,)li:«(aO 99.3l^:« 
2.2l°:°(me)l?:?(a0 SLStlH 


AS-^ 1 
Decay modes 


9k = 45° 
BRs (10-*^) LPFs (%) 


9k - -45° 
BRs (10-S) LPFs (%) 


Be ^ Ki(1270)"bi(1235)+ 
Be Ki(1400)%i(1235)+ 


^■G^oiimc)tli{a.) 91.31^:? 
1.3l°-^(TOe)l°-5(a») 100.0 ±0.0 


i.4l°:^(me)i;^:?(aO 100.01°:° 

1.5l°:«(m,)lJ:^(aO 93.61^:? 


Be ^ Ki(1270)+bi(1235)" 
Be ^ Ki(1400)+bi(1235)° 


0.8l°:^(m.)l°:°(a0 91.411? 
0.7l°:2(me)l°:^(aO 100.0 ±0.0 


o.7tVii^c)tVM) 100.01°:° 

0.8l°:^(me)l°:^(a0 93.6111 
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show us the dramatically different features: the former is dominated by the longitudi- 
nal components(~ 95%) while the latter governed by the transverse ones(~ 30%) when 
9k = 45°. The main reason is that the interferences induced by Be — >• K^j^cj) and 
Be K^^<t) are constructive(destructive) to B^ Ki{UOO)+ (piB^ Ki(1270)+(/)) in 
the two transverse polarizations, meanwhile, these interferences in the longitudinal po- 
larization contribute to these considered two decays oppositely. When 9k = —45°, the 
situation is quite the contrary. The decay mechanism and helicity structure for Be — ?■ K^cp 
decays will be tested by the LHCb and Super-B experiments. 

• For the AS = processes, Be — s- KiK^ modes, as presented in Table rvTIIl it is of interest 
to notice that the BRs for all these four considered decays are in the order of 10~^, 
which are within the reach of Be experiments at LHC greatly as discussed in Ref. [3]. 
These numerical results also present the strong dependance on the mixing angle 9k, 
which will also be tested by the relevant experiments in the near future. The longitudinal 
polarization fractions are around 95% ~ 100% within theoretical errors except for Be — s- 
;^(1400)°Ki(1400)+ (~ 73%) at 9k = 45° or Be K^imO)^ Ki{l270)+ (~ 72%) at 
9k = —45° and paly the dominant role. 

• As mentioned in the text above, although the suppressed CKM factor Vus ~ 0.22 is in- 
volved in the decay amplitudes (see Eqs. (|95I97I96I98P ) for the AS* = 1 Be ^ Ki{ai,bi) 
decays, the pQCD predictions of the BRs for Be — )• i^i^i are in the order of 10~^ and 
larger than that for Be — s- Kiai because the ^Pi meson behaves differently even contrarily 
to the ^Pi meson, whose behavior is close to that of the vector meson. For the polariza- 
tion fractions, all of these eight channels are governed by the longitudinal contributions 
evidently. From the numerical results as given in Table IVIIIt one can also find that the 
pQCD predictions of Be — )■ Kiai{Be — )• Kibi) are much(less) sensitive to the mixing angle 
9k- 

• For the AS = 1 Be ^ K+{fi{1285)Ji{U20)) decays, from the pQCD predictions 
presented in Table IIXI one can see that the contributions to the BRs for these four 
decays come from the overlap of various parts of Be — )■ K^^^fi, K^^fi, K^^fg, and 
-^lb/S) which have been given in Eqs. ()99til02p . Combining with four mixing parame- 
ters cos 9k, sm9K, cos^s and sin 6*3, these interferences result in the equivalent BRs for 



TABLE IX: Same as Table U but for 

6*3 = 38°(lst entry) and 6*3 = 50°(2nd entry). 



(i^i(1270)+,ii:i(1400)+)(/i(1285),/i(1420)) decays with 



AS-^ 1 
Decay modes 



Be ^ Ki(1270)+fi(1285) 
Be ^ Ki(1400)+fi(1285) 



Be ^ Ki(1270)+fi(1420) 
Be ^ Ki(1400)+fi(1420) 



45° 



BRs (10-^) 



LPFs (%) 



0.4("^c)_i 



1.7 



?(a, 
l{a, 



0.6 



+0.1^ 



Q 2(TOcy_0 



6 (a, 



4.3^i;:«(me)+;;?(a, 



4.4 



+0.5 



65.lt?^1 



69.1 



+22.1 

-19.6 



96.7t?i^5 

-4.0 
-12.8 



92.1; 



-fl6.9 



78.5 
71.9 



73.6±^| 



-1-4.8 
9.3 



-45° 



BRs (10- 



LPFs (%) 



1.6t°i(me)tl:J(a, 







+1 




lima): 



2 (a, 



-1-0.6 



4.4t° 
0.9t° 

-fO 





4(™c)ti: 



0.6 



5 

9(0, 

4/ 
7(0, 



4-2.7 

11.6 
-1-2.8 

-13.0 
4-27.2 
-19.4 



96.7 
92.1 
65.5 



71.5 



4-4.8 
-8.9 



70 94-4.8 
' '-'•^-9.3 
4-13.2 
34.4 
4-16.8 
-47.3 



81.9 
78.7 
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Be Ki(1270)+/i(1285) and Ki(1400)+/i(1285) decays, and suppressed one for 

Be iri(1270)+/i(1420) while enhanced one for Be i^i(1400)+/i(1420). Moreover, 
(a) the BRs for B^ R:^ fi{U20){Bc K^f i{1285)) depend strongly(weakly) on Ok 
for both 9s = 38° and 9^ = 50°; (b) the BRs for Be Ki(1270)+/i are more sensitive 
than that for Be — >■ -ftri(1400)^/i to 9^, when 9x = 45°, while the situation is quite the 
contrary when 9k = —45°; (c) the longitudinal contributions play an important role in 
all these considered channels. 

• Based on Eq. ([3]), apart from an overall sign, the physical states i<'i(1270) and 1^1(1400) 
can go one into another with changing the mixing angle 9k from 45° to —45° and vice 
versa, 

|i^l(1270))e^=45° = |i^l(1400)),^=_45°, 

|iri(1400))e^=45° = -\Ki{1270))e,,=-45°- (127) 

which further results in the decay amplitudes of Be — )■ Ki {V, A) (Here, A is a nonstrange 
axial- vector meson) as follows: 

A{Be^ Ki{mO)iV,A))e^=45^ =AiBe^Ki{im){V,A))e^=_45^, (128) 
A{Be^ Ki{im){V,A))e^=i5" = -^(S, ^ i^i(1270)(y, ^))e^=_45° • (129) 

These two relations, i.e., Eqs. (I128p and (I129p . can be manifested by the analytic formulas 
for Be Ki{V,A) decays as shown in Eqs. §MB, (ElEl]), I^MB and ([951 (106]) . The 
pQCD predictions for these considered Be — )• Ki{V,A) decays as listed in the second and 
third columns of Tables IVIIl IVIIH IIXI and |X] also display the phenomenologies induced 
by the same pattern. 

• For Be — )• KiK^ decays, however, it is not the shown in Eqs. (fT28]l and (fT29]) . 
According to the relation shown in Eq. (jl27p . there are some simple relations between the 
decay amplitudes as given in Eqs. (^VEM for Be KiKf decays : 

^(5,^;^(1270)i^i(1270))e^=45° = -^(5,^Kr(1400)Ki(1400))e^=_45° ,(130) 

^(Se^:^(1400)/Ci(1400))e^=45° = -^(^c ^ ^(1270)Ki(1270))e^=_45° ,(131) 

^(Se^;^(1270)i^i(1400))e^=45° = ^(S, ^ :^(1400)Ki(1270))e,,=_45° , (132) 

^(Se^:^(1400)/Ci(1270))e,,=45° = ^(S, ^ ;^(1270)Ki(1400))e^=_45° • (133) 



TABLE X: Same as Table U but for B, 
6*1 = 10° (1st entry) and Oi = 45° (2nd entry). 



(ii:i(1270)+,i^i(1400)+)(/ii(1170),/ii(1380)) decays with 



AS' = 1 
Decay modes 



Be ^ Ki(1270)+hi(1170) 
Be ^ Ki(1400)+hi(1170) 



Be ^ Ki(1270)+hi(1380) 
Be ^ Ki(1400)+hi(1380) 



45° 



BRs (10- 



LPFs (%) 



1.4 



+0.6 
-0.6 
^+0-3/ 



0.6t-^(me)t°i(a, 



-0.3*. 
^+0.7 



1.6^;]:^(me)+i;^(a, 
o(«' 



0.2 



0.5*. 
0.o("^c 



+0.2 



,+0.3 



0_9+o.3 



o.o(™f 

+0.5, 
-O.4I 



1.5+°-4 



i-0.3v"' 
\ + l. 



9(0, 



2.8tJi(m, 



'-0.9 

)_0.7('^' 



94.5 
87.9 



+2.3 
3.9 
+6.5 

-14.6 

93 n+^-^ 



98.5+°-« 

5+0.8 



-'-1.4 



98.6 
89 
94.4 



0.7 
+2.8 
-3.9 
+1.7 
-2.7 



-45° 



BRs (10- 



LPFs (%) 



i-6±g:I(m.): 





0.2;° 

i.4;« 

+0, 
-0 



0.5 



1. 



o(«.) 

o(a») 
lia^) 



98.5 
92.9 
94.6 
88.1 



+0.6 
0.9 
+7.5 

-15.1 
+2.3 
3.9 
+6.4 
-14.4 



1.5: 
2.8" 



-0.8 



0.9;° 



\+0. 



7(0*) 



[irir 



+2.9 
4.0 

+ 1.7 
-2.9 

98.51?:^ 
+0.9 

-0.5 



89.6 
94.3 



98.6 
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Of course, the above four relations, i.e., Eqs. (I130til33(l . can also be extracted from the 
pQCD predictions of BRs for Be — )• KiK^ decays as presented in Table rvTIIl apart from 
an overall sign. 

• At the first sight, it appears that the numerical results for Be — )• {fi,hi)(V, A)(Heve, 
A is either a or nonstrange axial-vector meson) decays are determined by the 
mixing angles ^3 and 9i, respectively, however, based on Ref. whose values will be 
eventually determined from 9k in Kia-Kib mixing system. Experimentally, it is thus 
very important to measure the channels precisely involving i('i(1270) and/or i^i(1400) to 
determine both of sign and size of the mixing angle Ok and reduce the uncertainties of 
theoretical predictions greatly. 

• The pQCD predictions for the CP-averaged branching ratios of considered Be decays vary 
in the range of 10~^ to 10^^. Since the LHC experiment can measure the Be decays 
with a branching ratio at 10~^ level 0], our pQCD predictions for the branching ratios 
of Be a+uj, hp, 'K^K+, 'KfK*+, /9+/i(1285), afb^, bfa^, a^/i(1285), a^/ii(1170), 
bfhi, K^K^ , ^t^i ^^"^ ^1^1 decays could be tested in the ongoing LHC experiments. 

• It is worth stressing that the theoretical predictions in the pQCD approach still have large 
theoretical errors induced by the still large uncertainties of many input parameters, e.g. 
Gegenbauer moments Oj. Any progress in reducing the error of input parameters, such as 
the Gegenbauer moments Oj and the charm quark mass nie, will help us to improve the 
precision of the pQCD predictions. 

V. SUMMARY 

In summary, we studied the sixty two charmless hadronic Be — )• VA, A A decays by employing 
the pQCD factorization approach based on the factorization theorem systematically. These 
considered decay channels can only occur via the annihilation type diagrams in the SM and 
they will provide an important platform for testing the magnitude and decay mechanism of the 
annihilation contributions and understanding the helicity structure of these considered channels 
and the content of the axial-vector mesons. Furthermore, these decay modes might also reveal 
the existence of exotic new physics scenario or nonperturbative QCD effects. 

The pQCD predictions for the CP-averaged branching ratios and longitudinal polarization 
fractions are displayed in Tables (jHIXl). From our perturbative evaluations and phenomenological 
analysis, we found the following results: 

• The pQCD predictions for the branching ratios vary in the range of 10^^ to 
10"^. There are many charmless Be — )• VA, AA decays with sizable branching 
ratios: Be afuj, bip, {Ki{1270)+ , Ki{im)+), (K^{1270f ,T<^iim)^)K*+ , 

p+/i(1285), a^lbl, b^al, a^/i(1285), a^hi{n70), 6+(/n(1170), /ii(1380)), 

(Kr(1270)°,;^(1400)°)(i^i(1270)+,i^i(1400)+), b^ {Ki{1260f , Ki{U00f) and 
(is:i(1270)+,i4:i(1400)+)(/ii(1170),/ii(1380)), which are with a decay rate at lO'^ 
or larger and could be measured at the LHC experiment. 

• For Be — )• VA, AA decays, the branching ratios of AS* = processes are generally much 
larger than those of AS* = 1 ones. Such differences are mainly induced by the CKM 
factors involved: Vud ~ 1 for the former decays while Vus ~ 0.22 for the latter ones. 
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• In general, since the behavior for ^Pi meson is much different from that for ^Pi meson, 
the branching ratios of the pure annihilation Be — )• Pi){V, Pi)) are larger than 
that of Be A{^Pi){V,A{^Pi)), which can be confronted with the LHC and Super-B 
experiments. 

• The longitudinal contributions play a dominant role in the most of these considered pure 
annihilation Be — )• VA, A A decays, which will be tested by the ongoing LHC and forth- 
coming Super-B experiments in the near future. 

• The pQCD predictions for several decays involving mixtures of ^Pi and/or ^Pi mesons 
are rather sensitive to the values of the mixing angles, which will be tested by the relevant 
experiments in the future. 

• Because only tree operators are involved, the CP-violating asymmetries for these consid- 
ered Be decays are absent naturally. 

• The pQCD predictions still have large theoretical uncertainties, mainly induced by the 
uncertainties of the Gegenbauer moments Oj in the meson distribution amplitudes. By 
reducing these uncertainties dramatically, one can improve the precision of the theoretical 
predictions effectively. 

• We here calculated the branching ratios and polarization fractions of the pure annihila- 
tion Be — )■ VA, AA decays by employing the pQCD approach. We do not consider the 
possible long-distance contributions, such as the rescattering effects, although they should 
be present, and they may be large and affect the theoretical predictions. It is beyond the 
scope of this work. 
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